Spectral properties of non-hydrogenic atoms in weak external fields 
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We study how the ionic core in a non-hydrogenic atom 
modifies the dynamics of a Rydberg electron in the presence 
of a weak static external field. We show that such a sys- 
tem is neither regular nor chaotic: its energy levels display 
unusual statistical properties, intermediate between the stan- 
dard Poisson and Random Matrix ones. The ionic core acts as 
a scatterer whose size is comparable to the de Broglie wave- 
length of the electron, inducing specific quantum effects. 

PACS: 05.45.+b, 32.60.+i, 03.65.Sq, 32.80.Rm 



The hydrogen atom in the presence of an external mag- 
netic field is a quantum system whose corresponding clas- 
sical dynamics is either regular or chaotic, depending on 
the scaled energy e = i?7~ 2 / 3 ' (E is the energy and 7 
the magnetic field in atomic units), and thus constitutes 
a prototype for studies on quantum chaos B| . 

In correspondance with the evolution of the classi- 
cal dynamics from quasi-integrable for a weak magnetic 
field (low negative scaled energy < —0.5) to chaotic at 
high field (high scaled energy > —0.13), the statistical 
properties of the energy levels evolve from a Poissonian 
spectrum to those of the Gaussian Orthogonal Ensem- 
ble (GOE) of random matrices jjj. In the first case, 
the nearest-neighbor spacing (NNS) distribution - which 
measures the probability to find the next energy level at 
distance s (in unit of the mean level spacing) - is P(s) = 
exp(— s), and in the second case, it is approximately given 
by the Wigner distribution P(s) = irs/2 exp(— its 2 /A). 

In a non-hydrogenic atom - for example an alkali atom 
as used in several experiments - the Rydberg elec- 
tron is in a highly excited state while all other electrons 
are in low excited states. The system (nucleus+inncr 
electrons) can be considered as a frozen ionic core with 
spherical symmetry, the size of few Bohr radii. When it 
is outside this ionic core, the Rydberg electron experi- 
ences a Coulomb field created by a charge Z=l at the 
nucleus. It is only when it penetrates the core that it 
feels a stronger force. At the scale of the Rydberg elec- 
tron (few thousands Bohr radii), the ionic core appears as 
a very small object perturbing the hydrogenic dynamics. 

A non-trivial question - addressed in this letter - is 
to determine the spectral properties of a non-hydrogenic 
atom in an external field. When the classical hydrogenic 
dynamics is chaotic, it seems that they remain described 
by the GOE pj. The situation is totally different when 
the hydrogenic dynamics is regular because the ionic core 



breaks the quasi-integrability. Courtney et al Q have 
shown a dramatic change in the statistical properties of 
energy levels which seem to be close to the GOE ones 
in the presence of the ionic core. They attribute this 
phenomenon to "core induced chaos" ; we here show that 
this point of view is incomplete. 

An efficient numerical method based on i?-matrix 
theory exists for finding the energy levels of a non- 
hydrogenic atom in an external field It relies on the 
physical picture outlined above: one splits space in a 
outer region where the effect of the ionic core is negligi- 
ble and where the Schrodinger equation for the Rydberg 
electron is solved by expansion on a suitable basis, and an 
inner region where the quantum dynamics is dominated 
by the ionic core (the effect of the external field is negli- 
gible) and can be described by a set of quantum defects. 
The matching between the solutions in the two regions 
gives the energy levels of the non-hydrogenic atom. 

We have calculated the lowest 40000 states of the hy- 
drogen atom in a magnetic field at scaled energy e = —0.5 
(for the L z = 0, even parity series) and the same set 
for the simplest non-hydrogenic atom, imposing a sin- 
gle non-zero quantum defect 5e=o = 0.5. This is not 
too far from a real atom like Lithium where 60 = 0.4 
and 5i>2 <C 0.01. Figure [I] shows the cumulative spacing 
distribution N(s) = P{x) dx obtained for various se- 
quences of energy levels (for o^ =0 = 0.5), all obtained at 
the same scaled energy, that is for the same hydrogenic 
dynamics. The lowest sequence is well described by the 
Wigner distribution, but, for higher excited states, the 
distribution deviates from Wigner and tends to a well 
defined limit which presents level repulsion - like the 
Wigner distribution - at the origin, but also has a much 
longer tail for large spacings - like the Poisson distri- 
bution, see Figs. |I|(d) and (e). Clearly, the behaviour 
at large s is exponential, not Gaussian. For comparison, 
we have also plotted the "semi-Poisson" distribution (de- 
fined below) fa] : 



P(s) = 4s exp(-2s), 



(1) 



which is in excellent agreement with our numerical re- 
sults. We have carefully checked that the numerically 
obtained distribution does not further change for higher 
excited states. In other words, when we go to higher and 
higher Rydberg states at fixed scaled energy, the statis- 
tical properties of the energy levels do not tend to the 
ones of the GOE. 
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FIG. 1. Statistical properties of the energy levels of a 
non-hydrogenic atom (quantum defect <5^ = o = 0.5) in a weak 
magnetic field (scaled energy -0.5, L z =0, even parity), (a) 
Cumulative nearest-neighbor spacing distribution between 
the 200 th and the 700" 1 excited states (solid histogram), 
compared with the Poisson (dashed line), semi-Poisson 
(dash-dotted line) and Wigner (dotted line) distributions; (b) 
and (c) Same as (a), respectively in the ranges 1000 th to 
2000 th and 5000 th to 40000 th excited states; (d) and (e) The 
nearest-neighbor spacing distribution itself, on linear and log- 
arithmic scales (data as in c). While, for low excited states, 
the distribution is well described by the Wigner distribution, 
the distribution for higher excited states is quite different, 
with a linear level repulsion like Wigner and an exponen- 
tial tail like Poisson. The agreement with the intermediate 
semi-Poisson distribution is excellent, see especially the tail 
at large spacing in (e) . In the corresponding hydrogenic situa- 
tion, all the distributions are close to the Poisson distribution. 

A similar phenomenon is observed in an external static 
electric field. For scaled energy e = EF^ 1 ! 2 (where F is 
the electric field in atomic units) sufficiently below the 
classical ionization threshold e = — 2, the highly excited 
states are quasi-discrete. In the absence of an ionic core, 
the system is integrable and displays a Poisson statistics. 
For a non-hydrogenic atom, the NNS distribution is lin- 
ear at small spacing (the cumulative distribution starts 
as s 2 ) and falls exponentially at large s - see Fig. |^. Al- 
though it slightly deviates from the semi-Poisson distri- 
bution (see discussion below), it has the same qualitative 
behaviour at small and large s and is clearly very different 
from a Wigner distribution. We expect such a behaviour 
to be general for non-hydrogenic atoms in any external 
field weak enough to keep the classical hydrogenic dy- 
namics quasi-integrable. 

A complete understanding requires a careful analysis 
of the physical phenomena taking place in the vicinity of 



the ionic core. The crucial point is to know whether a 
semiclassical approximation may be used there: we argue 
in the following that the answer is negative and that a 
specific approach is needed. 
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FIG. 2. Same as figure 1, but in a weak electric field (scaled 
energy -2.2) for states ranging from the 20000 th to the 40000 th 
excited states. Again, the cumulative spacing distribution is 
clearly different from the Poisson and the Wigner distribu- 
tions, in good qualitative agreement with the semi-Poisson 
prediction at small and large s. 

The ionic core acts as a spherical obstacle which scat- 
ters the Rydberg electron. The size of the scattering 
object (few Bohr radii) is comparable to the de Broglie 
wavelength of the Rydberg electron close to the core, 
even if one considers extremely highly excited states: the 
ionic core cannot be treated as a classical object || . The 
effect of a spherical scatterer is accounted for by a set 
of phase shifts in the different spherical channels. Here, 
these phase shifts are nothing but the quantum defects. 
Different model potentials for the ionic core may give 
the same set of quantum defects and consequently the 
same spectra in presence of an external field. They will 
however lead to different classical dynamics M and con- 
sequently to different thresholds for core-induced chaos: 
this is because near the nucleus, where the core strongly 
affects the classical trajectories and makes them to de- 
pend sensitively on the initial conditions, the semiclassi- 
cal approximation breaks down. It is thus impossible to 
build a solution of the Schrodinger equation which fol- 
lows closely one such classical trajectory, the wavepacket 
being unavoidably scattered near the origin ||. However, 
the breakdown of the semiclassical approximation is not 
extremely severe because the de Broglie wavelength is 
comparable - not much larger - to the size of the ionic 
core. A semiclassical approximation using core-scattered 
orbits thus give a correct qualitative understanding of 
the physics, for example the existence of additional mod- 
ulations in the density of states and density of oscillator 
strengths If quantitative results are needed - for ex- 
ample the amplitudes and phases of these modulations - 
a quantum ad hoc treatment of the ionic core is required, 
as done in B. From a temporal point of view, a typical 
classical trajectory follows an invariant torus of the hy- 
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drogenic motion, being from time to time scattered by 
the ionic core to another torus. The time interval be- 
tween two consecutive scattering events is of the order 
of the Heisenberg time, 7i/(mean level spacing), which 
expresses that the quantum dynamics can be considered 
neither as chaotic, nor as regular: it is an intermediate 
situation. 

The physics is very similar to the one of a rectangular 
two-dimensional billiard containing a small circular scat- 
terer, whose size is tuned to the de Broglie wavelength 
P). While, for a relatively large obstacle, the system (a 
standard Sinai billiard) has the statistical properties of 
the GOE, in the limiting case of a point scatterer, these 
are somewhat different: there is linear level repulsion at 
small spacing, but the tail significantly deviates from a 
Gaussian, exactly as our numerical observations on non- 
hydrogenic atoms. Because of particular attention on the 
level repulsion phenomenon, the distributions obtained in 
H have been incorrectly interpreted as Wigner distribu- 
tions and thus as proofs of the existence of "wave chaos" 
of "chaos induced by quantization" . 

The analogy between a non-hydrogenic atom and such 
a billiard can be made more formal. Suppose we are able 
to solve exactly the hydrogenic problem in the presence 
of an external field. Ei will denote the energy levels and 
■0i(r) the eigenstates. The Green's function is: 



(2) 



At any energy E, G(E, r,r') (viewed as a function of r) 
is a solution of the Schrodinger equation for the elec- 
tron, decreasing at infinity. Moreover, except in the im- 
mediate vicinity of the nucleus, it is also a solution of 
the Schrodinger equation for the non-hydrogenic atom. 
Hence, the energy spectrum of the non-hydrogenic atom 
in the presence of the external field is simply given by 
those values of E for which G(E, r, r') near r = can 
be matched with the field-free non-hydrogenic eigenfunc- 
tions (which depend on the quantum defects) described 
in |IoJ. In the specific case of a single non-zero quan- 
tum defect S — Se = o, only the £ = spherical component 
plays a role as the other ones cancel at the origin, and 
we finally end with the equation: 



E 



|^(0)| 2 = COt(TT^) 

Ei -E 2 ' 



(3) 



where the sum over i, which is formally divergent, has to 
be conveniently renormalized. This renormalization pro- 
cess is somewhat tricky and will be explained - together 
with the general case of several nonzero quantum defects 
- in detail elsewhere |ll]] . It is analogous to the one used 
in j^] for a point scatterer in a billiard (for which a simi- 
lar equation is obtained for the energy spectrum) , taking 
into account the presence of the Coulomb potential. 



Eq. (g) not only allows to calculate the non-hydrogenic 
energy levels as soon as the hydrogenic ones are known, 
but it also gives much insight in the statistical properties 
of the energy levels. A first remark is that there is ex- 
actly one root between two poles, i.e. one non-hydrogenic 
energy level in the interval limited by two consecutive hy- 
drogenic levels. We now concentrate on the case 6 = 0.5 
where the rhs in Eq. (|) is zero. The statistical properties 
of the roots of Eq. (B) can be obtained from the statisti- 
cal properties of the hydrogenic energy levels Ei and the 
numerators ^(O)! 2 . For a quasi-integrable situation, the 
Ei have a Poisson distribution, but the distribution of the 
numerators is not universal. It is then easy to show that 
the NNS distribution of the roots of Eq. (||) decreases 
like exp(— 2s) at large distance (unlike the Wigner distri- 
bution) and presents a linear level repulsion P(s) = as 
at small s with the coefficient a depending on the dis- 
tribution of the numerators. When all numerators are 
equal, one gets a = significantly smaller than 

the semi-Poisson prediction a — 4, see Eq. (0) ||. This 
is the situation realized for a weak external electric field 
where the eigenstates are separable in parabolic coordi- 
nates, and their overlaps with the usual spherical hydro- 
genic states are given by a Clebsch-Gordan coefficient, 
which is the same for all states in the case £ = m = : 
the slope near s — is smaller than 4, see Fig. ^, and 
a global small deviation from semi-Poisson is observed. 
For a non-uniform distribution of the numerators, a in- 
creases, being for usual smooth distributions close to 4, 
either smaller or larger. This is the case for an external 
magnetic field where the numerators are widely spread 
and the global agreement with semi-Poisson almost per- 
fect. In order to understand qualitatively the spectral 
properties, Bogomolny et al f| have introduced a sim- 
ple model - - baptized "semi-Poisson model" -for solving 
approximately Eq. (||), where each zero is lying exactly 
in the middle of two consecutive poles. The correspond- 
ing NNS distribution - called semi-Poisson distribution 
- is given by Eq. (0) and reproduces well our numerical 
results, see figs. |l| and |[ 

Several models based on different physical grounds pre- 
dict linear level repulsion. For example, the "short-range 
Dyson model" discussed in || (where there is level re- 
pulsion only between adjacent states, in contrast with 
the GOE where there is level repulsion between any pair 
of states) predicts a NNS distribution exactly equal to 
the semi-Poisson distribution, although other statistical 
quantities behave differently from the ones of the semi- 
Poisson model. In order to discriminate between the var- 
ious models, we study the next-nearest-neighbour distri- 
bution Pa(s) i.e. the statistical distribution of spacing 
between states n and n + 2. The corresponding cumula- 
tive distribution A^(s) = J P2{x) dx is shown in Fig. (j^) 
on a double logarithmic scale, in comparison with the pre- 
dictions of the Poisson, semi-Poisson, short-range Dyson 
and GOE models, which behave as s 2 , s 3 , s 4 and s 5 re- 
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spectively. Clearly, the numerical results scale as s 3 in 
good agreement with the semi-Poisson model, in sharp 
contrast with some other systems displaying intermedi- 
ate level statistics BH, which are closer to the short- 
range Dyson model. Note that semi-Poisson model does 
not predict correctly the coefficient of s 3 . Like the NNS, 
there is no universality here and the model is too crude. 

When the quantum defect S in the I = channel is not 
equal to 0.5, the rhs in Eq. (||) is not zero and, in average, 
the roots are not in the middle of consecutive poles. For 
equally spaced Ei and equal numerators, Eq. (0) can be 
exactly solved, and the roots lie exactly at fraction S 
of each interval. To describe the general situation (not 
equally spaced Ei), we define a 5— Poisson model where 
each level lies at SEi + (1 — S)E i+ i. The corresponding 
spacing distribution is: 




Spacing s 



know, non-hydrogenic atoms are the first "real" systems 
where this is numerically observed and where an experi- 
mental confirmation could be obtained in a near future. 
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for fruitful discussions. CPU time has been provided by 
IDRIS. Laboratoire Kastler-Brossel, de l'Ecole Normalc 
Superieure et de l'Universite Pierre et Marie Curie, is 
unite associee 18 du CNRS. 
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FIG. 4. Spacing distributions for a non-hydrogenic atom 
in a weak magnetic field (scaled energy -0.5, from 13000* to 
22000 th excited states) with quantum defect 8e=o = 0.15 (dot- 
ted line) and 0.85 (dashed line) in comparison with the analyt- 
ical prediction of the 5-Poisson model, Eq. (^). The excellent 
agreement validates the model which smoothly evolves from a 
Poisson distribution for 8 = to a semi-Poisson distribution 
at 8 — 0.5 where the effect of the ionic core is maximum. 



FIG. 3. Cumulative next-nearest-neighbor spacing distri- 
butions (on a double logarithmic scale) for a non-hydrogenic 
atom in a weak magnetic (data of figure lc) or electric 
(data of figure 2) field, compared with the Poisson (dashed 
line), semi-Poisson (solid line), short-range Dyson model 
(dash-dotted line) and GOE (dotted line) distributions. The 
functional dependance at small spacing s is clearly propor- 
tional to s 3 , in agreement with the semi-Poisson model, and 
deviates from all other distributions. 

It is shown in Fig. (ji|) in comparison with the NNS 
distributions obtained numerically for a non-hydrogenic 
atom in a magnetic field with 5 = 0.15 and 0.85. The 
agreement is excellent. The same distribution is obtained 
for quantum defects 5 and 1 — S, as predicted by Eq. (Ji|). 

In conclusion, we have shown that, for a non- 
hydrogenic atom in a weak external field, the presence 
of the ionic core which scatters the Rydbcrg electron, 
leads to statistical properties of the energy levels of a 
new type, intermediate between the standard Poisson 
and GOE statistics. Similar observations have been re- 
cently reported in other systems |||l2],|l3| , which indi- 
cates their importance and broad interest: as far as we 
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